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ABSTRACT.  present  a  summary  of  results  concerning  the  determina¬ 
tion  of  an  unknown  conductivity  by  means  of  static  measurements  at 
the  boundary.  The  main  emphasis  is  on  identif lability ;  fie^only 
briefly  discuss  the  reconstruction  problem.  Some  references  are  given 
to  related  work  for  time  dependent  problems. 


1.  INTRODUCTION.  We  study  the  following  inverse  problem:  can  one  determine 
an  unknown  conductivity  y(x)  Inside  a  body  ft  ,  by  means  of  static  measurements 
at  the  boundary?  A.  P.  Calderon  raised  this  question  in  [2];  it  may 
be  seen  as  a  natural  extension  of  a  problem  analyzed  by  Cannon,  Douglas  and 
Jones  in  1963  [3,4].  Despite  some  progress, many  unsolved  problems  remain. 

Our  main  goal  here  is  to  summarize  what  is  known  about  identifiabllity;  we 
do  this  in  sections  2  and  3,  which  are  based  mostly  on  [18].  A  few  of  the 
results  presented  -  notably  2E  and  3C  -  are  previously  unpublished.  In 
section  4,  we  touch  on  the  reconstruction  problem,  given  finitely  many  measure¬ 
ments;  and  section  5  reviews  the  literature  on  some  related  problems.  Poten¬ 
tial  applications  include  nondestructive  testing  and  water  resources  management 
[14,13],  but  these  will  not  be  discussed  here. 

Throughout,  ft  will  be  a  bounded  domain  in  ®n  ,  with  unit  normal  v  along 
3ft  .  The  unknown  conductivity  y(x)  may  take  real  scalar  or  matrix  values, 
corresponding  to  an  Isotropic  or  anisotropic  material.  In  the  isotropic  case 

(1.1)  Y  €  L*(ft)  ,  ess  inf  y  >  0  , 

while  In  the  anisotropic  case 

(1.2)  YtJ  -  Yji  €  L*(ft)  ,  XU|2  «  (y(*)C,5)  S  AU|2 

for  all  x  f  ft  and  C  €  *n  ,  with  A  >  0  .  We  consider  solutions  of 

(1.3)  Iyj  ■  V*(y(x)Vu)  ■  0  ; 

in  the  context  of  heat  conduction  u  represents  temperature  and  yVu  the  heat 

*  ~ 
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flux;  in  the  context  of  direct  current  electrical  conduction  u  represents 
voltage  and  y7u  is  the  vector  of  current  flow.  The  natural  things  to  measure 

at  3ft  are  the  Dirichlet  data  uL_  and  the  Neumann  data  y(x)7u*v|  .  We 

1/2  -1/2  ““  ° 
denote  by  P  :  H  (3ft)  -*•  H  '  (3ft)  the  operator  which  associates  the  former 

to  the  latter, 


we  shall  say  that 


Knowledge  of 
(1.4) 


PY<t>  -  yVu- v  1 3Q  with  L^u  -  0  ,  u|3f}  -  *  ; 

'\1  and  ^2  give  the  same  boundary  measurements"  if 


P^  yields  the  energy  quadratic  form 


qy<*> 


(y7u, 

Jn 


7u)dx 


4i*P  ids  , 
3ft  Y 


by  Green's  formula,  where  the  rightmost  integral  represents  the  dual  pairing 
1/2  -1/2 

of  H  and  H  .  Conversely,  Q  determines  P  by  the  polarization 

identity;  hence  "y.  and  y-  give  the  same  boundary  measurements”  iff 
1  z  1/2 

Qv  (41)  ■  Qv  (<t>)  for  each  €  H  (3ft)  .  The  following  variational  charac- 

T1  t2 

ter izat ions  of  are  well-known: 

(1.5)  Qv($)  -  min  [  (yVw,7w)dx 

w  €  HX(ft)  ift 
w"$  on  3ft 


(1.6) 


-Qy<*) 


min 

<xa2(«;Rn 

V*o«0 


>( 


(y-1o,a)dx 

ft 


$a-vds 


2.  IDENTIFIABILITY  -  THE  ISOTROPIC  CASE.  In  one  dimension,  only  the  harmonic 
mean  of  y  can  be  detected  by  boundary  measurement:  it  is  an  easy  exercise 
to  show  that 

2A.  For  ft  ■  (a,b)  ,  y^  and  Y2  give  the  same  boundary  measurements  iff 
they  have  the  same  harmonic  mean. 

Fortunately,  the  situation  is  entirely  different  in  dimension  greater  than  one. 

Cannon,  Douglas,  and  Jones  considered  cylindrical  domains,  with  y  con¬ 
stant  along  lines  parallel  to  the  axis,  in  1963.  They  showed  that  such  y 
are  identifiable: 

2B  [3].  Suppose  ft  ■  G  x  (a,b)  ,  with  G  a  bounded.  C^,flt  domain  in  ®n  ^  , 
a  >  0  .  Then  the  elements  of 

rx  -  (y  €  C1,a(G)  :  inf  y  >  0} 

can  be  distinguished  by  means  of  boundary  measurements. 

Their  procedure  for  reconstructing  y  is  remarkably  direct.  Taking 

(a,b)  ■  (0,ir)  for  simplicity,  and  writing  x'  •  (x^ . . .  .x„_1)  ,  let  L^u  ■  0 
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u»0  on  G*{0,ir} 

u  -  g(x')sin  xn  on  3G  *  (0,ir)  , 

where  g  :  30  -»■  H  is  any  positive  function.  If 

k(x')  -  y(x')  £-\ 

3xn 1 xn-0 

is  measured,  then 


y  -  k-exp(-w)  , 


where  w(x’)  solves 


k  on  G 
tin  g  on  3G 


Thus  y(~')  is  determined  for  all  x'  ,  using  a  single  choice  of  g  . 

The  restriction  that  y  be  independent  of  x^  is,  of  course,  crucial  to 
the  preceeding  analysis.  One  is  not  really  finding  y  "in  the  interior", 
since  it  is  determined  by  its  values  along  the  boundary.  When  the  dependence 
of  y  is  unrestricted,  one  naturally  obtains  information  at  the  boundary  more 
easily  than  in  the  interior.  If  everything  is  smooth,  then  y  is  determined 
to  infinite  order  at  30  . 

2C  [18].  Suppose  that  30  is  smooth,  and  that  y1»Y2  e  c  (0)  •  If.  Yj_  and 
Y2  Klve  the  same  boundary  measurements,  then 

(2.1)  -  D— y2  on  30 

for  all  k  -  (k^, . . .  ,kR)  >,  0  ,  where 


The  proof  of  2C  is  local  in  character,  but  not  constructive.  For  xQ  €  30 

with  v  (xn)  0  ,  consider  the  Dirichlet  data 
n  U 

y  -  1  n-1 

♦N(x)  -IT  n  0  (N(Xj-XQ  .j)) 

with  corresponding  solutions  u^  , 

v*-°-  ■  *- •  ‘“•2- 

If  *  «  Cj  OR)  has  vanishing  moments  of  order  $  M-l  ,  then  a  version  of 
"St.  Vsnant's  principle"  provides  that 

|7uJ(x)|  (  C  N-M  ,  N  -  « 

for  x  f  0  bounded  away  from  xQ  . 

If  (2.1)  fails,  then  (relabeling  if  necessary) 
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Y^x)  -  y2(x)  *  c  diet  (x,3ft)^ 

with  ^  ;  0  and  C  >  0  ,  in  a  neighborhood  of  some  Xq  €  3ft  ;  it  now  follows 
that  *  ^Y2^N^  provided  M  >  nt/2  and  N  is  sufficiently  large,  a 

contradiction  to  the  fact  that  y^  and  y2  give  the  same  boundary  measure¬ 
ments.  Details  may  be  found  in  [18];  see  also  3C  below. 

An  immediate  corollary  of  2C  is  this: 

2D  :  For  <i  smoothly  bounded  domain  ft  c  Rn  ,  n  z  2  ,  the  elements  of 

r2  »  {restrictions  to  ft  of  positive,  real 

analytic  functions  defined  in  a  neighborhood 
of  ft  } 


can  be  distinguished  by  boundary  measurements. 


The  analogue  of  2D  with  "analytic"  replaced  by  "piecewise  analytic"  is 
open.  We  are  encouraged,  however,  by  the  following  example. 

2 

2E.  Let  ft  be  the  unit  disc  in  R  ,  with  polar  coordinates  (r,0)  ,  and 
denote  by  the  set  of  conductivities 


y(r,9)  -  |Y°  r  ^ 

l*  1  rn  i  r  <  1 


with  0  <  Tq  <  1  and  y^  a  positive  constant .  Then  the  elements  of  can 

be  distinguished  by  boundary  measurements . 


To  prove  2E  ,  consider  y,y  f  ,  with  y  corresponding  to  yQ  ,  fQ  ;  we 
shall  show  that  Q^(sin  NO)  Q~(sin  N8)  for  all  sufficiently  large  N  ,  unless 
y  «  y  .  If  fg  «  Tq  this  follows  instantly,  since  either  y  <  y  or  y  >  y 
and  consequently  either  Q-(sin  N6)  <  Q^Csin  N0)  or  Q*(sln  N0)  >  Q^(sin  N0)  , 
provided  yQ  +  YQ  • 

Relabeling  if  necessary,  we  assume  that  fg  <  rQ  .  For  N  %  0  ,  let  u^ 
solve 

Ly  “n  "  °  ’  "  8in  N0  ’ 


and  notice  that  for  r  <  rg 

(2-2>  “u  -  cNrNsin  **  ■  CN  -  iW7.  ~~~T  • 

r0  (l-y0)+(y0+l) 


In  case 

(2.3) 


Yq  >  1  ,  (2.2)  implies  that 

.  Y0l7uN|2dx  <  L  (Yg-l)l%|2dx 
ro  ro<r<ro 


for  sufficiently  large  N,  whence 
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*dx  -  [  Yq  I  I  2dx  +  [  |VuN|2dx  + 

Jrcrn  JTn<r<rn  Jr 


|VuN|Mx 


0  0 


'ro<r<ro 


71^1 2dx  +  [  I^Ujjl^x  , 

Jr„<r<l 


|  yjVUjj^dx  <  |  ylVu^dx  -  Q^,(sin  N0) 


It  follows,  using  (1.5),  that  Q~(sin  N0)  <  Q^(sin  N0)  . 

Next,  suppose  that  Yq  <  1  •  In  this  case,  a  similar  argument  gives 


f  Y~1iaijl2dx  <  f  y-1IcnI 

Jn  <*52 


with  oN  ■  yVujj  ,  and  N  sufficiently  large.  Using  the  dual  variational 
principle  (1.6),  we  conclude  that  Q-(sin  N0)  >  Q^(sin  N0)  . 

3.  IDENTIFIABILITY  -  THE  ANISOTROPIC  CASE.  In  the  anisotropic  case,  one  can 
not  expect  to  recover  the  full  matrix  y^  ,  as  the  following  two  examples 
demonstrate. 


3A.  Let  52  C  R  ,  n  *  1  ,  and  let  y 


(1.2).  For  any  C  dlffeomor- 


phism  *  :  0  -*■  SI  with 


*(x)  *  x  ,  D4(x)  -  I  for  all  x  £  352  , 


y^(4(x))  -  |  det(D$ (x))|  i*D4(x)t  *y(x)  *Di»(x)  . 


Then  all  elements  of 


r4  ■  (y*  :  $  satisfies  (3.1)} 


live  the  same  boundary  measurements. 


Me  owe  this  remark  to  L.  Tartar.  If  L  u  -  0  ,  then  L  .  u  *  0  ,  with 

Y  (y  ) 

u*(x)  *  u®$  l(x)  ; 

by  (3.1),  u*  ■  u  on  y*»Vu*«  y*7u  on  352  . 

2 

3B  [25].  Let  52  be  the  unit  disc  in  R  ,  with  polar  coordinates  (r,0)  . 
For  any  function  a(r)  ,  let 

(a  cos20  +  a  ^sin20  (a-a  ^)sin  0*cos  0  \ 

-1  2-12J' 

(a-a  )sin  9*cos  0  a  sin  0  +  a  cos  9  / 


Then  mil  elements  of 


Tj  ■  (y°  :  a  <  L  (0,1)  ,  ess  inf  a  >  0} 
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give  the  same  boundary  measurements ■ 

Indeed , 

r-L  a  -  (ra(r)  T~)+  — Vr- 
ya  3r  \  3r  /  ra(r) 


39 


note  that  when  a  =  1  this  is  just  r  times  the  Laplacian.  For  N  €  Z  ,  the 
solution  of 

IN  6 


has  the  form  v(r)e 


This  implies  that 


iN9 


V  “N  ■  0  ’  “n| 


sn 


with  v(l)  »  1  and 


'h  (ra(r) 


3v\  _  N 
3r  /  ra(r) 


0  . 


v(r)  -  cx  exp  (iNlj^  +  c2  exp  (- 1 H | jgfc)  , 

with  +  c_  *  1  .  Since  v(r)eiN9  €  H^Q)  ,  c 2  must  equal  zero;  hence  the 


Neumann  data  associated  to 


“N 


is 


r-1 


3uN  ]  |  IN 6 

“IT"  |N'e 


1.N0  1/2 

regardless  of  the  choice  of  a  .  The  span  of  (e  }  is  dense  in  H  '  (3fl)  , 
so  each  ya  £  T^  gives  the  same  boundary  measurements. 

What  can  one  detect,  in  the  anisotropic  case?  The  natural  analogue  of  2C 
is  this:  if  (n-1)  eigenvalues  and  eigenvectors  of  y  are  known,  the  last 
eigenvalue  can  be  distinguished  by  boundary  measurements. 


3C.  Let  y,y  be  two  symmetric ,  positive  definite  matrices  with  entries  in 
L  (fl)  ,  and  let  { A^> , { J  and  te^Jjte^}  be  the  corresponding  eigenvalues 
and  eigenvectors.  For  Xq  €  3fl  ,  let  B  be  ii  neighborhood  of  Xq  relative 


to  U  , 

and 

suppose  that 

(3.2) 

Y.Y 

€  C*(B) 

,  and  3fJ  fl  B 

la 

00 

C  ; 

(3.3) 

Sj  "  6j 

’  V 

Xj  in  B  ,  1  : 

?  J  * 

n-1 

(3.4) 

e 

n 

(xQ)  *v  (Xq)  +  0 

H1/2 

(3.5) 

QyW  - 

Q-(0 

for  every  <t>  € 

(an) 

with 

supp  4> 

c  B  0  38  . 

Then 

(3.6) 

n(V  "  D^^n(x0) 

for  any 

k  - 

0*1  *  •  •  •  ^ 

0  . 
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We  sketch  the  proof.  For  a  fixed  z  €  90  near  Xg  ,  let  a 

sequence  of  functions  on  90  such  that 

(3.7) 

(3.8) 


supp  d>N  +  { Z  } 


Mk.annB  *  ckN  *for  a11  k5'M 


(3.9) 


90ns 


i  , 


where  the  norms  are  standard  Sobolev  norms  based  on  L  ;  the  existence  of 
such  a  sequence  may  be  deduced  from  [18],  for  any  fixed  M  >  0.  If  solves 

S'V  ’  0  in  n  *  “n  "  0n  ’ 

then  (3. 7)-(3.9)  imply 

-M 


(3.10) 


5  CN 


for  all  neighborhoods  U  of  z  (c.f.  Lemma  2  of  [18]).  Condition  (3.4)  gives 
in  a  sufficiently  small  neighborhood  U  of  z 


(3.11) 


’  2  2 
| w |  ds  $  C  |e  *Vw|  dx 

ungo  'u  n 


for  any  w  €  H  (ft)  with  w  *  0  in  U  ;  using  this  and  (3.7)-(3.9),  one 
obtains 


Z, 


ir-(n+l+e)Z 


•  ' 0 


0.12)  P  IvVuJ  **  * 

for  any  Z  >  0  ,  e  >  0  ,  with  p(x)  “  dist  (x,90)  ,  provided  M  >  (n+l)£/2 
(cf.  Lemma  3  of  [18]). 

In  order  to  prove  (3.6)  it  is  sufficient  to  verify  that  Y  *  (^r)  Y 

in  a  90  -  neighborhood  of  Xq  ,  for  any  k  >,  0  .  We  prove  this  by  contradic¬ 
tion,  using  (3.10)  and  (3.12).  If  it  fails,  we  may  assume  (switching  y  and 
y  if  necessary)  that  there  exists  a  z  near  Xq  and  a  neighborhood  U  C  B  of 
z  such  Chat 

(3.13) 

with  C  >  0  ,  Z  >.  0  .  Then 


\  (x)  -  \  (x)  >A  Cp(x)  for  x  ?  U  , 
n  n 


!  5  (,  rt,v7Vd* 

+  C  JD°£[’=-7uNl2d-  ■ 

using  (3.3)  and  (3.13).  If  M  >  (n+l)Z/2  ,  then  (3.10)  and  (3.12)  show  that 

C  V7uN|2dx  ”  lnVu(97U»’7Vdx 

for  large  K.  Hence,  using  (1.5), 

W  “  ln(Y7uN*7VdX  ”  |n(*7V7Vdx  *  Qy(V  ; 

This  contradicts  (3.5),  and  the  proof  is  complete. 
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The  preceding  argument  used  (3.4)  to  justify  (3.11).  The  following 
example  suggests  that  (3.4)  may  be  dispensable. 


3D  [4].  For  0  -  (a,b)  x  (c,d)  c  J  ,  let  denote  the  family  of  conduc- 

tivlties 


Y(x) 


a(xzV 


for  a  €  C  with  inf  a  >  0  .  Then  the  elements  of  can  be  distinguished 

by  boundary  measurements.  More  specifically,  if  y  €  ,  i  =  1,2  and 


1  2 

y^7u  ‘m  -  y^7u  * v  at  x^  -  b 


where  u  is  the  solution  to 

i 


Ly_u  “0  in  0 


0  at  x. 


a,b 


u  “  sin 


(■=6) 


at  x_ 


c,d 


then 


Y1  =  Y2  • 


4.  RECONSTRUCTION.  Setting  aside  the  question  of  identif lability,  how  might 


one  estimate  y  in  practice?  A  straightforward  approach  is  the  following: 


let  1/  be  a  finite-dimensional  subspace  of  H^O)  j  for  w^  6  V  , 

1/2  3ul 

f  H  (30)  ,  and  measure  y  >  where  u  solves 

30  3v  j 


1  {  j  s  m  ,  set  -  Wj 


L  u.  »  0  in  0  ,  u 
Y  j  j 


30 


If  G  is  a  finite-parameter  family  of  possible  conductivities,  and  y  €  G  , 


let  iij  €  V  be  the  solution  to  the  Galerkin  equation 


[  (yVu  ,V$)dx  ■  0  for  all  $  €  V  f!  S^(0) 

Jo  3 


u  •  6^  on  30 


-  3w  -i/2 

Now  assume  furthermore  that  V  and  G  are  selected  so  that  y  ^  *  (30) 


whenever  w  €  V  and  y  €  G  ,  and  choose  y  to  minimise 

3u 


(4.1) 


J(r)  -  l  IIy 
J-l 


3v  Y  3v 


3u,  2 

^  H“1/2(30) 


among  all  y  €  C  • 
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A  variation  of  this  method,  and  its  finite-difference  analogue,  is  studied 
by  Falk  in  (11),  for  operators  of  the  form 

-  yu"  +  c(x)u  =*  f  (x)  0  <  x  <  1 

c,f  known;  y  >  0  constant 

in  one  dimension.  He  takes  G  »  H+  ,  and  uses  just  one  measurement  (m*l)  at 
one  point  on  the  boundary  (x»0),  so  that  the  functional  *  0  for  some  y  ;  and 
he  estimates  how  |y-y|  depends  on  the  choice  of  !/  . 

The  minimization  of  (4.1)  has  apparently  not  been  studied  in  higher 
dimensions.  There  is,  however,  a  large  literature  on  parameter  identification, 
much  of  which  is  relevant:  see,  for  example,  (10,19,22).  In  other  problems 
where  identif icability  and  well-posedness  are  known,  one  can  often  prove  the 
convergence  of  such  a  procedure  (1).  And  in  some  cases,  even  identifiability 
can  be  proved  using  an  approximation  algorithm  (20). 

Calderdn  takes  a  completely  different  approach  to  the  reconstruction 
problem  in  [2],  for  the  scalar  case  with  <S  »  !|y-1||  sufficiently  small.  He 
shows  that  the  Fourier  transform  of  y  (extended  by  zero  off  ft)  has  the  form 

Y<S)  -  f(5)  +  R(S)  , 

where  f(£)  can  be  determined  by  boundary  measurements,  and 

|R(£)  I  s  05^  >exp(ff  - 1  £  |  •diam(ft)  )  . 

Thus  boundary  measurements  suffice  to  approximate  y(S)  well  at  low  frequencies, 
if  5  is  small.  For  any  fixed  £  €  Rn  ,  measuring  f(£)  requires  the  Neumann 
data  from  just  one  solution  of  L^u  «  0  ,  corresponding  to  Dirichlet  conditions 
exp[ir(i£+n)  *x)  ,  with  n  €  ,  n*5  -  0  ,  |n)  *  1 5 i  . 

5.  RELATED  WORK.  The  parabolic  analogue  of  our  problem  is  to  determine 
Y  *  y(*>t)  >  given  overdetermined  boundary  data  for  solutions  of 

|jr  -  V*(yVu)  -  0 

on  a  space-time  cylinder  ft  x  (0,T)  .  This  problem  h.-'S  been  studied  extensive¬ 
ly  in  space-dimension  one,  both  as  to  the  identifiability  of  y  and  as  to  its 
numerical  approximation:  see  [6]  for  y  »  constant,  [5,9,15,16]  for  y  *  Y(t)  , 
and  [17,21,26]  for  y  -  y(x)  •  We  know  of  no  results  in  space  dimension 
greater  than  one  for  spatially  varying  y  . 

An  interesting  nonlinear  analogue  is  obtained  by  letting  y  depend  on  u  , 
so  that  the  equation  becomes  7*(y(x,u)Vu)  -  0  .  The  case  y  -  y(u)  has  been 
studied  in  [7],  and  a  related  parabolic  problem  is  treated  in  [8]. 

Many  authors  have  studied  the  reconstruction  of  an  unknown  y(x)  .  given 
knowledge  of  a  single  function  u  everywhere  on  ft  ,  satisfying  V'(yVu)  •  0  . 
This  problem  is  of  particular  interest  for  studying  ground-water  flow  through 
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porous  rock.  Ident if lability  is  analyzed  in  [24],  and  the  convergence  of 
numerical  schemes  are  studied  in  [23,12].  Applications  and  other  numerical 
methods  have  been  discussed  in  a  dozen  or  so  articles  in  Water  Resources 
Research  over  the  last  ten  years,  of  which  [13,27]  are  examples. 
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